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SIGNAL PROCESSING WITH CHANNEL EIGENMODE 
DECOMPOSITION AND CHANNEL INVERSION FOR MIMO 
SYSTEMS 

BACKGROUND 

Field 

[1001] The present invention relates generally to data communication, and more 
specifically to techniques for performing signal processing with channel eigenmode 
decomposition and channel inversion for multiple-input multiple-output (MIMO) 
communication systems. 

Background 

[1002] A multiple-input multiple-output (MIMO) communication system employs 
multiple (Nt) transmit antennas and multiple (N R ) receive antennas for data 
transmission. A MTMO channel formed by the Nt transmit and Nr receive antennas 
may be decomposed into Ns independent channels, with N s < min {N T ,N R }. Each of 
the Ns independent channels is also referred to as a spatial subchannel of the MIMO 
channel and corresponds to a dimension. The MIMO system can provide improved 
performance (e.g., increased transmission capacity) if the additional dimensionalities 
created by the multiple transmit and receive antennas are utilized. 
[1003] The spatial subchannels of a wideband MIMO system may encounter 
different channel conditions due to various factors such as fading and multipath. Each 
spatial subchannel may thus experience frequency selective fading, which is 
characterized by different channel gains at different frequencies (i.e., different 
frequency bins or subbands) of the overall system bandwidth. With frequency selective 
fading, each spatial subchannel may achieve different signal-to-noise-and-interference 
ratios (SNRs) for different frequency bins. Consequently, the number of information 
bits per modulation symbol (or data rate) that may be transmitted at different frequency 
bins of each spatial subchannel for a particular level of performance (e.g., 1% packet 
error rate) may be different from bin to bin. Moreover, because the channel conditions 
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typically vary with time, the supported data rates for the bins of the spatial subchannels 
also vary with time. 

[1004] To combat frequency selective fading in a wideband channel, orthogonal 
frequency division multiplexing (OFDM) may be used to effectively partition the 
system bandwidth into a number of (N F ) subbands (which may also be referred to as 
frequency bins or subchannels). With OFDM, each frequency subchannel is associated 
with a respective subcarrier upon which data may be modulated. For a MTMO system 
that utilizes OFDM (i.e., a MIMO-OFDM system), each frequency subchannel of each 
spatial subchannel may be viewed as an independent transmission channel. 
[1005] A key challenge in a coded communication system is the selection of the 
appropriate data rates and coding and modulation schemes to be used for a data 
transmission based on channel conditions. The goal of this selection process is to 
maximize throughput while meeting quality objectives, which may be quantified by a 
particular packet error rate (PER), certain latency criteria, and so on. 
[1006] One straightforward technique for selecting data rates and coding and 
modulation schemes is to "bit load" each transmission channel in the MIMO-OFDM 
system according to its transmission capability, which may be quantified by the 
channel's short-term average SNR. However, this technique has several major 
drawbacks. First, coding and modulating individually for each transmission channel 
can significantly increase the complexity of the processing at both the transmitter and 
receiver. Second, coding individually for each transmission channel may greatly 
increase coding and decoding delay. And third, a high feedback rate would be needed 
to send channel state information (CSI) indicative of the channel conditions (e.g., the 
gain, phase, and SNR) of each transmission channel. 

[1007] For a MIMO system, transmit power is another parameter that may be 
manipulated to maximize throughput. In general, the overall throughput of the MIMO 
system may be increased my allocating more transmit power to transmission channels 
with greater transmission capabilities. However, allocating different amounts of 
transmit power to different frequency bins of a given spatial subchannel tends to 
exaggerate the frequency selective nature of the spatial subchannel. It is well known 
that frequency selective fading causes inter-symbol interference (ISI), which is a 
phenomenon whereby each symbol in a received signal acts as distortion to subsequent 
symbols in the received signal. The ISI distortion degrades performance by impacting 
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the ability to correctly detect the received symbols. To mitigate the deleterious effects 
of ISI, equalization of the received symbols would need to be performed at the receiver. 
Thus, a major drawback in frequency-domain power allocation is the additional 
complexity at the receiver to combat the resultant additional ISI distortion. 
[1008] There is therefore a need in the art for techniques to achieve high overall 
throughput in a MfMO system without having to individually code each transmission 
channel and which mitigate the deleterious effects of ISI. 

SUMMARY 

[1009] Techniques are provided herein for processing a data transmission at a 
transmitter and a receiver of a MIMO system such that high performance (e.g., high 
overall throughput) is achieved. In an aspect, a time-domain implementation is 
provided which uses frequency-domain channel eigen-decomposition, channel 
inversion, and (optionally) "water-pouring" results to derive pulse-shaping and beam- 
steering solutions for the transmitter and receiver. 

[1010] Channel eigen-decomposition is performed at the transmitter to determine 
the eigenmodes (i.e., the spatial subchannels) of a MEMO channel and to obtain a first 
set of steering vectors, which are used to precondition modulation symbols prior to 
transmission over the MIMO channel. Channel eigen-decomposition may be performed 
based on an estimated channel response matrix, which is an estimate of the (time- 
domain or frequency-domain) channel response of the MIMO channel. Channel eigen- 
decomposition is also performed at the receiver to obtain a second set of steering 
vectors, which are used to condition received symbols such that orthogonal symbol 
streams are recovered at the receiver. _ 
[1011] Channel inversion is performed at the transmitter to derive weights, which 
are used to minimize or reduce the amount of ISI distortion at the receiver. In 
particular, the channel inversion may be performed for each eigenmode used for data 
transmission. One set of weights may be derived for each eigenmode based on the 
estimated channel response matrix for the MEMO channel and is used to invert the 
frequency response of the eigenmode. 

[1012] Water-pouring analysis may (optionally) be used to more optimally allocate 
the total available transmit power to the eigenmodes of the MEMO channel. In 
particular, eigenmodes with greater transmission capabilities may be allocated more 
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transmit power, and eigenmodes with transmission capabilities below a particular 
threshold may be omitted from use (e.g., by allocating these bad eigenmodes with zero 
transmit power). The transmit power allocated to each eigenmode then determines the 
data rate and possibly the coding and modulation scheme to be used for the eigenmode. 
[1013] At the transmitter, data is initially processed (e.g., coded and modulated) in 
accordance with a particular processing scheme to provide a number of modulation 
symbol streams (e.g., one modulation symbol stream for each eigenmode). An 
estimated channel response matrix for the MIMO channel is obtained (e.g., from the 
receiver) and decomposed (e.g., in the frequency domain, using channel eigen- 
decomposition) to obtain a set of matrices of right eigen-vectors and a set of matrices of 
singular values. A number of sets of weights are then derived based on the matrices of 
singular values, with each set of weights being used to invert the frequency response of 
a respective eigenmode used for data transmission. Water-pouring analysis may also be 
performed based on the matrices of singular values to obtain a set of scaling values 
indicative of the transmit powers allocated to the eigenmodes. A pulse-shaping matrix 
for the transmitter is then derived based on the matrices of right eigen-vectors, the 
weights, and the scaling values (if available). The pulse-shaping matrix comprises 
steering vectors, which are used to precondition the streams of modulation symbols to 
obtain streams of preconditioned symbols to be transmitted over the MIMO channel. 
[1014] At the receiver, the estimated channel response matrix is also obtained (e.g., 
based on pilot symbols sent from the transmitter) and decomposed to obtain a set of 
matrices of left eigen-vectors. A pulse-shaping matrix for the receiver is then derived 
based on the matrices of left eigen-vectors and used to condition a number of received 
symbol streams to obtain a number of recovered symbol streams. The recovered 
symbols are further processed (e.g., demodulated and decoded) to recover the 
transmitted data. 

[1015] Various aspects and embodiments of the invention are described in further 
detail below. The invention further provides methods, digital signal processors, 
transmitter and receiver units, and other apparatuses and elements that implement 
various aspects, embodiments, and features of the invention, as described in further 
detail below. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
[1016] The features, nature, and advantages of the present invention will become 
more apparent from the detailed description set forth below when taken in conjunction 
with the drawings in which like reference characters identify correspondingly 
throughout and wherein: 

[1017] FIG. 1 is a block diagram of an embodiment of a transmitter system and a 
receiver system in a MIMO system; 

[1018] FIG. 2 is a block diagram of a transmitter unit within the transmitter system; 
[1019] FIGS. 3A and 3B are diagrams that graphically illustrate the derivation of 
the weights used to invert the frequency response of each eigenmode of a MIMO 
channel; 

[1020] FIG. 4 is a flow diagram of a process for allocating the total available 
transmit power to the eigenmodes of the MIMO channel; 

[1021] FIGS. 5 A and 5B are diagrams that graphically illustrate the allocation of the 
total transmit power to three eigenmodes in an example MIMO system; 
[1022] FIG. 6 is a flow diagram of an embodiment of the signal processing at the 
transmitter unit; 

[1023] FIG. 7 is a block diagram of a receiver unit within the receiver system; and 
[1024] FIG. 8 is a flow diagram of an embodiment of the signal processing at the 
receiver unit. 

DETAILED DESCRIPTION 
[1025] The techniques described herein for processing a data transmission at a 
transmitter and receiver may be used for various wireless communication systems. For 
clarity, various aspects and embodiments of the invention are described specifically for 
a multiple-input multiple-output (MIMO) communication system. 
[1026] A MIMO system employs multiple (N T ) transmit antennas and multiple (N R ) 
receive antennas for data transmission. A MIMO channel formed by the Nt transmit 
and Nr receive antennas may be decomposed into Ns independent channels, with 
N s < min {N T , N R } . Each of the N s independent channels is also referred to as a spatial 
subchannel of the MIMO channel. The number of spatial subchannels is determined by 
the number of eigenmodes for the MIMO channel, which in turn is dependent on a 
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channel response matrix that describes the response between the N T transmit and N R 
receive antennas. 

[1027] FIG. 1 is a block diagram of an embodiment of a transmitter system 1 10 and 
a receiver system 150, which are capable of implementing various signal processing 
techniques described herein. 

[1028] At transmitter system 1 10, traffic data is provided from a data source 1 12 to 
a transmit (TX) data processor 114, which formats, codes, and interleaves the traffic 
data based on one or more coding schemes to provide coded data. The coded traffic 
data may then be multiplexed with pilot data using, for example, time division multiplex 
(TDM) or code division multiplex (CDM), in all or a subset of the data streams to be 
transmitted. The pilot data is typically a known data pattern processed in a known 
manner, if at all. The multiplexed pilot and coded traffic data is interleaved and then 
modulated (i.e., symbol mapped) based on one or more modulation schemes to provide 
modulation symbols. In an embodiment, TX data processor 114 provides one 
modulation symbol stream for each spatial subchannel used for data transmission. The 
data rate, coding, interleaving, and modulation for each modulation symbol stream may 
be determined by controls provided by a controller 130. 

[1029] The modulation symbols are then provided to a TX MIMO processor 120 
and further processed. In a specific embodiment, the processing by TX MIMO 
processor 120 includes (1) determining an estimated channel frequency response matrix 
for the MIMO channel, (2) decomposing this matrix to determine the eigenmodes of the 
MIMO channel and to derive a set of "steering" vectors for the transmitter, one vector 
for the modulation symbol stream to be transmitted on each spatial subchannel, (3) 
deriving a transmit spatio-temporal pulse-shaping matrix based on the steering vectors 
and a weighting matrix indicative of the transmit powers assigned to the frequency bins 
of the eigenmodes, and (4) preconditioning the modulation symbols with the pulse- 
shaping matrix to provide preconditioned modulation symbols. The processing by TX 
MIMO processor 120 is described in further detail below. Up to N T streams of 
preconditioned symbols are then provided to transmitters (TMTR) 122a through 122t. 
[1030] Each transmitter 122 converts a respective preconditioned symbol stream 
into one or more analog signals and further conditions (e.g., amplifies, filters, and 
frequency upcon verts) the analog signals to generate a modulated signal suitable for 
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transmission over the MIMO channel. The modulated signal from each transmitter 122 
is then transmitted via a respective antenna 124 to the receiver system. 
[1031] At receiver system 150, the transmitted modulated signals are received by N R 
antennas 152a through 152r, and the received signal from each antenna 152 is provided 
to a respective receiver (RCVR) 154. Each receiver 154 conditions (e.g., filters, 
amplifies, and frequency downconverts) the received signal, digitizes the conditioned 
signal to provide a stream of samples, and further processes the sample stream to 
provide a stream of received symbols. An RX MTMO processor 160 then receives and 
processes the N R received symbol streams to provide Nt streams of recovered symbols, 
which are estimates of the modulation symbols transmitted from the transmitter system. 
In an embodiment, the processing by RX MIMO processor 160 may include (1) 
determining the estimated channel frequency response matrix for the MIMO channel, 
(2) decomposing this matrix to derive a set of steering vectors for the receiver, (3) 
deriving a receive spatio-temporal pulse-shaping matrix based on the steering vectors, 
and (4) conditioning the received symbols with the pulse-shaping matrix to provide the 
recovered symbols. The processing by RX MIMO processor 160 is described in further 
detail below. 

[1032] A receive (RX) data processor 162 then demodulates, deinterleaves, and 
decodes the recovered symbols to provide decoded data, which is an estimate of the 
transmitted traffic data. The processing by RX MTMO processor 160 and RX data 
processor 162 is complementary to that performed by TX MIMO processor 120 and TX 
data processor 114, respectively, at transmitter system 1 10. 

[1033] RX MIMO processor 160 may further derive channel impulse responses for 
the MIMO channel, received noise power and/or signal-to-noise-and-interference ratios 
(SNRs) for the spatial subchannels, and so on. RX MTMO processor 160 would then 
provide these quantities to a controller 170. RX data processor 162 may also provide 
the status of each received packet or frame, one or more other performance metrics 
indicative of the decoded results, and possibly other information. Controller 170 then 
derives channel state information (CSI), which may comprise all or some of the 
information received from RX MTMO processor 160 and RX data processor 162. The 
CSI is processed by a TX data processor 178, modulated by a modulator 180, 
conditioned by transmitters 154a through 154r, and sent back to transmitter system 110. 
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[1034] At transmitter system 110, the modulated signals from receiver system 150 
are received by antennas 124, conditioned by receivers 122, and demodulated by a 
demodulator 140 to recover the CSI transmitted by the receiver system. The CSI is then 
provided to controller 130 and used to generate various controls for TX data processor 
1 14 and TX MIMO processor 120. 

[1035] Controllers 130 and 170 direct the operation at the transmitter and receiver 
systems, respectively. Memories 132 and 172 provide storage for program codes and 
data used by controllers 130 and 170, respectively. 

[1036] Techniques are provided herein for achieving high performance (e.g., high 
overall system throughput) via a time-domain implementation that uses frequency- 
domain channel eigen-decomposition, channel inversion, and (optionally) water-pouring 
results to derive time-domain pulse-shaping and beam-steering solutions for the 
transmitter and receiver. 

[1037] Channel eigen-decomposition is performed at the transmitter to determine 
the eigenmodes of the MEMO channel and to derive a first set of steering vectors, which 
are used to precondition the modulation symbols. Channel eigen-decomposition is also 
performed at the receiver to derive a second set of steering vectors, which are used to 
condition the received symbols such that orthogonal symbol streams are recovered at 
the receiver. The preconditioning at the transmitter and the conditioning at the receiver 
orthogonalize the symbol streams transmitted over the MIMO channel. 
[1038] Channel inversion is performed at the transmitter to flatten the frequency 
response of each eigenmode (or spatial subchannel) used for data transmission. As 
noted above, frequency selective fading causes intersymbol interference (ISI), which 
can degrade performance by impacting the ability to correctly detect the received 
symbols at the receiver. Conventionally, the frequency selective fading may be 
compensated for at the receiver by performing equalization on the received symbol 
streams. However, equalization increases the complexity of the receiver processing. 
With the inventive techniques, the channel inversion is performed at the transmitter to 
account for the frequency selective fading and to mitigate the need for equalization at 
the receiver. 

[1039] Water-pouring (or water-filling) analysis is used to more optimally allocate 
the total available transmit power in the MIMO system to the eigenmodes such that high 
performance is achieved. The transmit power allocated to each eigenmode may then 
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determine the data rate and the coding and modulation scheme to be used for the 
eigenmode. 

[1040] These various processing techniques are described in further detail below. 
[1041] The techniques described herein provide several potential advantages. First, 
with time-domain eigenmode decomposition, the maximum number of eigenmodes with 
different SNRs is given by min (N T ,N R ) . If one independent data stream is transmitted 
on each eigenmode and each data stream is independently processed, then the maximum 
number of different coding/modulation schemes is also given by min(N T , N R ) . It is 
also possible to make the received SNRs for the data streams essentially the same, 
thereby further simplifying the coding/modulation. The techniques described herein can 
thus greatly simplify the coding/modulation for a data transmission by avoiding the per- 
bin bit allocation required to approach channel capacity in MIMO-OFDM systems that 
utilize frequency-domain water-pouring. 

[1042] Second, the channel inversion at the transmitter results in recovered symbol 
streams at the receiver that do not require equalization. This then reduces the 
complexity of the receiver processing. In contrast, other wide-band time-domain 
techniques typically require complicated space-time equalization to recover the symbol 
streams. 

[1043] Third, the time-domain signaling techniques described herein can more 
easily integrate the channel/pilot structures of various CDMA standards, which are also 
based on time-domain signaling. Implementation of the channel/pilot structures may be 
more complicated in OFDM-based systems that perform frequency-domain signaling. 
[1044] FIG. 2 is a block diagram of an embodiment of a transmitter unit 200, which 
is capable of implementing various processing techniques described herein. Transmitter 
unit 200 is an embodiment of the transmitter portion of transmitter system 110 in FIG. 
1. Transmitter unit 200 includes (1) a TX data processor 114a that receives and 
processes traffic and pilot data to provide Nt modulation symbol streams and (2) a TX 
MTMO processor 120a that preconditions the modulation symbol streams to provide NT- 
preconditioned symbol streams. TX data processor 1 14a and TX MIMO processor 120a 
are one embodiment of TX data processor 114 and TX MIMO processor 120, 
respectively, in FIG. 1. 
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[1045] In the specific embodiment shown in FIG. 2, TX data processor 114a 
includes an encoder 212, a channel interleaver 214, and a symbol mapping element 216. 
Encoder 212 receives and codes the traffic data (i.e., the information bits, d,) in 
accordance with one or more coding schemes to provide coded bits. The coding 
increases the reliability of the data transmission. In an embodiment, a separate coding 
scheme may be used for the information bits for each eigenmode (or spatial subchannel) 
selected for use for data transmission. In alternative embodiments, a separate coding 
scheme may be used for each subset of spatial subchannels, or a common coding 
scheme may be used for all spatial subchannels. The coding scheme(s) to be used are 
determined by controls from controller 130 and may be selected based on the CSI 
received from the receiver system. Each selected coding scheme may include any 
combination of cyclic redundancy check (CRC), convolutional coding, Turbo coding, 
block coding, and other coding, or no coding at all. 

[1046] Channel interleaver 214 interleaves the coded bits based on one or more 
interleaving schemes. Typically, each selected coding scheme is associated with a 
corresponding interleaving scheme. The interleaving provides time diversity for the 
coded bits, permits the data to be transmitted based on an average SNR of each spatial 
subchannel used for the data transmission, combats fading, and further removes 
correlation between coded bits used to form each modulation symbol. 
[1047] Symbol mapping element 216 then receives and multiplexes pilot data with 
the interleaved data and further maps the multiplexed data in accordance with one or 
more modulation schemes to provide modulation symbols. A separate modulation 
scheme may be used for each spatial subchannel selected for use, or for each subset of 
spatial subchannels. Alternatively, a common modulation scheme may be used for all 
selected spatial subchannels. 

[1048] The symbol mapping for each spatial subchannel may be achieved by 
grouping sets of bits to form data symbols (each of which may be a non-binary value) 
and mapping each data symbol to a point in a signal constellation corresponding to the 
modulation scheme selected for use for that spatial subchannel. The selected 
modulation scheme may be QPSK, M-PSK, M-QAM, or some other scheme. Each 
mapped signal point is a complex value and corresponds to a modulation symbol. 
Symbol mapping element 216 provides a vector of modulation symbols for each symbol 
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period, with the number of modulation symbols in each vector corresponding to the 
number of spatial subchannels selected for use for that symbol period. Symbol mapping 
element 216 thus provides up to Nt modulation symbol streams. These streams 
collectively form a sequence of vectors, with are also referred to as the modulation 
symbol vectors, s(n) , with each such vector including up to N s modulation symbols to 
be transmitted on up to Ns spatial subchannels for the n-th symbol period. 
[1049] Within TX MTJvlO processor 120a, the response of the MTMO channel is 
estimated and used to precondition the modulation symbols prior to transmission to the 
receiver system. In a frequency division duplexed (FDD) system, the downlink and 
uplink are allocated different frequency bands, and the channel responses for the 
downlink and uplink may not be correlated to a sufficient degree. For the FDD system, 
the channel response may be estimated at the receiver and sent back to the transmitter. 
In a time division duplexed (TDD) system, the downlink and uplink share the same 
frequency band in a time division multiplexed manner, and a high degree of correlation 
may exist between the downlink and uplink channel responses. For the TDD system, 
the transmitter system may estimate the uplink channel response (e.g., based on the pilot 
transmitted by the receiver system on the uplink) and may then derive the downlink 
channel response by accounting for any differences such as those between the transmit 
and receive antenna array manifolds. 

[1050] In an embodiment, the channel response estimates are provided to TX 
MEVIO processor 120a as a sequence of N R xN T matrices, A(n) , of time-domain 
samples. This sequence of matrices is collectively referred to as a channel impulse 
response matrix, A . The (t, j) -th element, A, ; , of the estimated channel impulse 
response matrix, A, for i = (l, 2, ... ,N R ) and ; = (1, 2, ... ,N T ), is a sequence of 
samples that represents the sampled impulse response of the propagation path from the 
j-th transmit antenna to the i-th receive antenna. 

[1051] Within TX MEVIO processor 120a, a fast Fourier transformer 222 receives 
the estimated channel impulse response matrix, A (e.g., from the receiver system) and 
derives the corresponding estimated channel frequency response matrix, H, by 
performing a fast Fourier transform (FFT) on the matrix (i.e., H = FFT [fi] ). This 
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may be achieved by performing an Nf- point FFT on a sequence of N F samples for each 
element of A to derive a set of N F coefficients for the corresponding element of H , 
where Nf corresponds to the number of frequency bins for the FFT (i.e., the length of 
the FFT). The N R ■ N T elements of H are thus N R • N T sets of coefficients representing 
the frequency responses of the propagation paths between the Nt transmit antennas and 
N R receive antennas. Each element of H is the FFT of the corresponding element of 
A. The estimated channel frequency response matrix, H, may also be viewed as 
comprising a set of N F matrices, R(k) for k = (0, 1, ... , N F — 1) . 

Channel Eigen-Decomposition 
[1052] A unit 224 then performs eigen-decomposition of the MIMO channel used 
for data transmission. In one embodiment for performing channel eigen-decomposition, 
unit 224 computes the singular value decomposition (SVD) of the estimated channel 
frequency response matrix, H . In an embodiment, the singular value decomposition is 
performed for each matrix U(k), for fc = (0, 1, ... ,N F -l). The singular value 
decomposition of matrix H(fc) for frequency bin k (or frequency f k ) may be expressed 
as: 

E(k) = V(k)A(k)Y."(k) , Eq(l) 

where U(ifc) is an N R x N R unitary matrix (i.e., U H U = I, where I is the identity 
matrix with ones along the diagonal and zeros everywhere else); 
A(k) is an N R xN T diagonal matrix of singular values of H(&) ; and 
V(fc) is an N T xN T unitary matrix. 

The diagonal matrix A(&) contains non-negative real values along the diagonal (i.e., 
A(^) = diag (A, (A:), A 2 (k), ... ,A„ T (k))) and zeros elsewhere. The A t (k), for 
i = (l, 2, ... ,N T ), are referred to as the singular values of the matrix H(&). The 
singular value decomposition is a matrix operation known in the art and described in 
various references. One such reference is a book by Gilbert Strang entitled "Linear 
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Algebra and Its Applications," Second Edition, Academic Press, 1980, which is 
incorporated herein by reference. 

[1053] The result of the singular value decomposition is three sets of N F matrices, 
U, A, and V" , where U = [U(0) ... U(£) ... U(JV F -1)] , and so on. For each value of 
it, LP) is the N R xN R unitary matrix of left eigen-vectors of H(fc), V(fc) is the 
N T xN T unitary matrix of right eigen-vectors of H(&), and A(fc) is the N R xN T 
diagonal matrix of singular values of H(£) . 

[1054] In another embodiment for performing channel eigen-decomposition, unit 
224 first obtains a square matrix R(k) as R(k) = u" (k)H(k) . The eigenvalues of the 
square matrix R(k) would then be the squares of the singular values of the matrix 
H(fc), and the eigen-vectors of R(k) would be the right eigen-vectors of H(fc), or 
\(k). The decomposition of R(k) to obtain the eigenvalues and eigen-vectors is 

known in the art and not described herein. Similarly, another square matrix R (k) may 
be obtained as R (k) = H(k)U H (k) . The eigenvalues of this square matrix R (k) 
would also be the squares of the singular values of H(&) , and the eigen-vectors of 
R (k) would be the left eigen-vectors of H(£) , or U(&) . 

[1055] The channel eigen-decomposition is used to decompose the MIMO channel 
into its eigenmodes, at frequency f k , for each value of k where k = (0, 1, ... , N F -1) . 
The rank r(k) of H(&) corresponds to the number of eigenmodes for the MIMO 
channel at frequency f k , which corresponds to the number of independent channels 
(i.e., the number of spatial subchannels) available in frequency bin k. 
[1056] As described in further detail below, the columns of Y(k) are the steering 
vectors associated with frequency f k to be used at the transmitter for the elements of 
the modulation symbol vectors, s(«) . Correspondingly, the columns of U(&) are the 
steering vectors associated with frequency f k to be used at the receiver for the elements 
of the received symbol vectors, r(«) . The matrices U(&) and V(fc) , for 
k = (0, 1, ... , N F -1) , are used to orthogonalize the symbol streams transmitted on the 
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eigenmodes at each frequency f k . When these matrices are used to precondition the 
modulation symbol streams at the transmitter and to condition the received symbol 
streams at the receiver, either in the frequency domain or the time domain, the result is 
the overall orthogonalization of the symbol streams. This then allows for separate 
coding/modulation per eigenmode (as opposed to per bin), which can greatly simplify 
the processing at the transmitter and receiver. 

[1057] The elements along the diagonal of Mk) are A u (k), for 
i = {l, 2, .... r(k)}, where r(k) is the rank of K(k) . The columns of V(k) and \(k), 
u ( (it) and Xi(k), respectively, are solutions to the eigen equation, which may be 
expressed as: 

i(k)y l (k) = JL a u t (k) . Eq (2) 

[1058] The three sets of matrices, U(£), A(Jfc), and V(*) , for 
k = (0, 1, ... , N F - 1) , may be provided in two forms - a "sorted" form and a "random- 
ordered" form. In the sorted form, the diagonal elements of each matrix A(k) are 
sorted in decreasing order so that ^(k) > A 22 (k) > ... > X rr (k) , and their eigen-vectors 
are arranged in corresponding order in V(k) and V(£) . The sorted form is indicated 
by the subscript s, i.e., U s (it) , A, (k) , and V s (k) , for k = (0, 1, ... , N F - 1) . 
[1059] In the random-ordered form, the ordering of the singular values and eigen- 
vectors may be random and further independent of frequency. The random form is 
indicated by the subscript r. The particular form selected for use, sorted or random- 
ordered, influences the selection of the eigenmodes for use for data transmission and the 
coding and modulation scheme to be used for each selected eigenmode. 
[1060] A weight computation unit 230 receives the set of diagonal matrices, A, 
which contains a set of singular values (i.e., ^{k), A 22 (k), ...,A rr (k)) for each 
frequency bin. Weight computation unit 230 then derives a set of weighting matrices, 
W, where W = [W(0) ... W(fc) ... W(N F -1)]. The weighting matrices are used to 
scale the modulation symbol vectors, s(n) , in either the time or frequency domain, as 
described below. 
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[1061] Weight computation unit 230 includes a channel inversion unit 232 and 
(optionally) a water-pouring analysis unit 234. Channel inversion unit 232 derives a set 
of weights, w . , for each eigenmode, which is used to combat the frequency selective 
fading on the eigenmode. Water-pouring analysis unit 234 derives a set of scaling 
values, b, for the eigenmodes of the MIMO channel. These scaling values are 
indicative of the transmit powers allocated to the eigenmodes. Channel inversion and 
water-pouring are described in further detail below. 

Channel Inversion 

[1062] FIG. 3A is a diagram that graphically illustrates the derivation of the set of 
weights, w l7 , used to invert the frequency response of each eigenmode. The set of 
diagonal matrices, A(k) for k = (0, 1, ... , N F -1) , is shown arranged in order along an 
axis 310 that represents the frequency dimension. The singular values, A ti (k) for 
i = (1, 2, ... ,N S ), of each matrix A(k) are located along the diagonal of the matrix. 
Axis 312 may thus be viewed as representing the spatial dimension. Each eigenmode of 
the MIMO channel is associated with a set of elements, {/t u (k)} for 
k = (0, 1, ... , N F -1) , that is indicative of the frequency response of that eigenmode. 
The set of elements {A. u (k)} for each eigenmode is shown by the shaded boxes along a 
dashed line 314. For each eigenmode that experiences frequency selective fading, the 
elements {A u (k) } for the eigenmode may be different for different values of k. 
[1063] Since frequency selective fading causes ISI, the deleterious effects of ISI 
may be mitigated by "inverting" each eigenmode such that it appears flat in frequency at 
the receiver. The channel inversion may be achieved by deriving a set of weights, 
{w a (k)} for k = (0, 1, ... ,N F -1), for each eigenmode such that the product of the 
weights and the corresponding eigenvalues (i.e., the squares of the diagonal elements) 
are approximately constant for all values of k, which may be expressed as 
w u (it) • X\ (k) = a, , for k = (0, 1, ... , N F - 1) . 

[1064] For eigenmode i, the set of weights for the N F frequency bins, 
w .. = [w.. (0) ... w.. (k) ... w„. (N F - l)f , used to invert the channel may be derived as: 
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w, (k) = , for k = (0, 1, ... , N F - 1) , 



Eq(3) 



where a,, is a normalization factor that may be expressed as: 




Eq(4) 



As shown in equation (4), a normalization factor a, is determined for each eigenmode 
based on the set of eigenvalues (i.e., the squared singular values), {Al(k)} for 
k = (0, 1, ... ,N F -l), associated with that eigenmode. The normalization factor a, is 

N F -l N F -l 

defined such that £ w u (k) = £ ^- (*) • 



[1065] FIG. 3B is a diagram that graphically illustrates the relationship between the 
set of weights for a given eigenmode and the set of eigenvalues for that eigenmode. For 
eigenmode i, the weight w u (k) for each frequency bin is inversely related to the 
eigenvalue Al(k) for that bin, as shown in equation (3). To flatten the spatial 
subchannel and minimize or reduce ISI, it is undesirable to selectively eliminate 
transmit power on any frequency bin. The set of N F weights for each eigenmode is used 
to scale the modulation symbols, s(n), in the frequency or time domain, prior to 
transmission on the eigenmode. 

[1066] For the sorted order form, the singular values A u (k) , for i = (1, 2, ... ,N S ), 
for each matrix A(k) are sorted such that the diagonal elements of A(k) with smaller 
indices are generally larger. Eigenmode 0 (which is often referred to as the principle 
eigenmode) would then be associated with the largest singular value in each of the iv> 
diagonal matrices, A(k) , eigenmode 1 would then be associated with the second largest 
singular value in each of the Nf diagonal matrices, and so on. Thus, even though the 
channel inversion is performed over all Nf frequency bins for each eigenmode, the 
eigenmodes with lower indices are not likely to have too many bad bins (if any). Thus, 
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at least for eigenmodes with lower indices, excessive transmit power is not used for bad 
bins. 

[1067] The channel inversion may be performed in various manners to invert the 
MEVIO channel, and this is within the scope of the invention. In one embodiment, the 
channel inversion is performed for each eigenmode selected for use. In another 
embodiment, the channel inversion may be performed for some eigenmodes but not 
others. For example, the channel inversion may be performed for each eigenmode 
determined to induce excessive ISI. The channel inversion may also be dynamically 
performed for some or all eigenmodes selected for use, for example, when the MIMO 
channel is determined to be frequency selective (e.g., based on some defined criteria). 
[1068] Channel inversion is described in further detail in U.S. Patent Application 
Serial No. 09/860,274, filed May 17, 2001, U.S. Patent Application Serial No. 
09/881,610, filed June 14, 2001, and U.S. Patent Application Serial No. 09/892,379, 
filed June 26, 2001, all three entitled "Method and Apparatus for Processing Data for 
Transmission in a Multi-Channel Communication System Using Selective Channel 
Inversion," assigned to the assignee of the present application and incorporated herein 
by reference. 

Water-Pouring 

[1069] In an embodiment, water-pouring analysis is performed (if at all) across the 
spatial dimension such that more transmit power is allocated to eigenmodes with better 
transmission capabilities. The water-pouring power allocation is analogous to pouring a 
fixed amount of water into a vessel with an irregular bottom, where each eigenmode 
corresponds to a point on the bottom of the vessel, and the elevation of the bottom at 
any given point corresponds to the inverse of the SNR associated with that eigenmode. 
A low elevation thus corresponds to a high SNR and, conversely, a high elevation 
corresponds to a low SNR. The total available transmit power, P tota u is then "poured" 
into the vessel such that the lower points in the vessel (i.e., those with higher SNRs) are 
filled first, and the higher points (i.e., those with lower SNRs) are filled later. A 
constant P set is indicative of the water surface level for the vessel after all of the total 
available transmit power has been poured. This constant may be estimated initially 
based on various system parameters. The power allocation is dependent on the total 
available transmit power and the depth of the vessel over the bottom surface. The 
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points with elevations above the water surface level are not filled (i.e., eigenmodes with 
SNRs below a particular value are not used for data transmission). 
[1070] In an embodiment, the water-pouring is not performed across the frequency 
dimension because this tends to exaggerate the frequency selectivity of the eigenmodes 
created by the channel eigenmode decomposition described above. The water-pouring 
may be performed such that all eigenmodes are used for data transmission, or only a 
subset of the eigenmodes is used (with bad eigenmodes being discarded). It can be 
shown that water-pouring across the eigenmodes, when used in conjunction with the 
channel inversion with the singular values sorted in descending order, can provide near- 
optimum performance while mitigating the need for equalization at the receiver. 
[1071] The water-pouring may be performed by water-pouring analysis unit 234 as 
follows. Initially, the total power in each eigenmode is determined as: 



where a 2 is the received noise variance, which may also be denoted as the received 
noise power N 0 . The received noise power corresponds to the noise power on the 
recovered symbols at the receiver, and is a parameter that may be provided by the 
receiver to the transmitter as part of the reported CSI. 

[1073] The transmit power, />, to be allocated to each eigenmode may then be 
determined as: 




Eq(5) 



[1072] The SNR for each eigenmode may then be determined as: 



Pi A 

SNR,. =-4 , 

/T 



Eq (6) 




Eq(7a) 




Eq(7b) 
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where P sel is a constant that may be derived from various system parameters, and P lotal 

is the total transmit power available for allocation to the eigenmodes. 

[1074] As shown in equation (7a), each eigenmode of sufficient quality is allocated 

transmit power of | P set — 1 . Thus, eigenmodes that achieve better SNRs are 



allocated more transmit powers. The constant P sel determines the amounts of transmit 
power to allocate to the better eigenmodes. This then indirectly determines which 
eigenmodes get selected for use since the total available transmit power is limited and 
the power allocation is constrained by equation (7b). 

[1075] Water-pouring analysis unit 234 thus receives the set of diagonal matrices, 
A, and the received noise power, a 2 . The matrices A are then used in conjunction 
with the received noise power to derive a vector of scaling values, 
b = [b 0 ... b t ... b N$ ] 7 , where b t = /> for i = (1, 2, ... ,N S ). The P, are the solutions to 
the water-pouring equations (7a) and (7b). The scaling values in b are indicative of the 
transmit powers allocated to the Ns eigenmodes, where zero or more eigenmodes may 
be allocated no transmit power. 

[1076] FIG. 4 is a flow diagram of an embodiment of a process 400 for allocating 
the total available transmit power to a set of eigenmodes. Process 400, which is one 
specific water-pouring implementation, determines the transmit powers, P t , for i G / , to 
be allocated to the eigenmodes in set /, given the total transmit power, P ma u available at 
the transmitter, the set of eigenmode total powers, P i X , and the received noise power, 



[1077] Initially, the variable n used to denote the iteration number is set to one (i.e., 
n = l) (step 412). For the first iteration, set I{n) is defined to include all of the 
eigenmodes for the MIMO channel, or /(n) = {1, 2, ... , N s } (step 414). The 
cardinality (or length) of set I(n) for the current iteration n is then determined as 
L, (n) = \l(n)\ , which is L, (n) = N s for the first iteration (step 416). 
[1078] The total effective power, P eff (n) , to be distributed across the eigenmodes in 
set I(n) is next determined (step 418). The total effective power is defined to be equal 
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to the total available transmit power, P lotal , plus the sum of the inverse SNRs for the 
eigenmodes in set I(n) . This may be expressed as: 

M») = C+E-^ • Eq(8) 

[1079] The total available transmit power is then allocated to the eigenmodes in set 
I(n). The index i used to iterate through the eigenmodes in set I(n) is initialized to 
one (i.e., i = 1) (step 420). The amount of transmit power to allocate to eigenmode i is 
then determined (step 422) based on the following: 

Each eigenmode in set I(n) is allocated transmit power, P. , in step 422. Steps 424 and 
426 are part of a loop to allocate transmit power to each of the eigenmodes in set I(n) . 
[1080] FIG. 5A graphically illustrates the total effective power, P eg , for an 
example MMO system with three eigenmodes. Each eigenmode has an inverse SNR 
equal to a 2 1 A\ , for i = {1, 2, 3} , which assumes a normalized transmit power of 1.0. 
The total transmit power available at the transmitter is P wsal = P l +P 2 +P 3 , and is 
represented by the shaded area in FIG. 5A. The total effective power is represented by 
the area in the shaded and unshaded regions in FIG. 5A. 

[1081] For water-pouring, although the bottom of the vessel has an irregular 
surface, the water level at the top remains constant across the vessel. Likewise and as 
shown in FIG. 5A, after the total available transmit power, P tota u has been distributed to 
the eigenmodes, the final power level is constant across all eigenmodes. This final 
power level is determined by dividing P eJJ (n) by the number of eigenmodes in set I(n) , 
L, (n) . The amount of power allocated to eigenmode i is then determined by 
subtracting the inverse SNR of that eigenmode, <J 2 /A*, from the final power level, 
P eff (n)/L, (ri) , as given by equation (9) and shown in FIG. 5A. 

[1082] FIG. 5B shows a situation whereby the water-pouring power allocation 
results in an eigenmode receiving negative power. This occurs when the inverse SNR 



WO 2004/002047 



PCT7US2003/019464 



21 

of the eigenmode is above the final power level, which is expressed by the condition 
(P eff (n)/L,(n))<(a 2 /Al). 

[1083] Referring back to FIG. 4, at the end of the power allocation, a determination 
is made whether or not any eigenmode has been allocated negative power (i.e., P t ; < 0) 
(step 428). If the answer is yes, then the process continues by removing from set I(n) 
all eigenmodes that have been allocated negative powers (step 430). The index n is 
incremented by one (i.e., n = n + 1 ) (step 432). The process then returns to step 416 to 
allocate the total available transmit power among the remaining eigenmodes in set I(n) . 
The process continues until all eigenmodes in set 7(n) have been allocated positive 
transmit powers, as determined in step 428. The eigenmodes not in set I(n) are 
allocated zero power. 

[1084] Water-pouring is also described by Robert G. Gallager, in "Information 
Theory and Reliable Communication," John Wiley and Sons, 1968, which is 
incorporated herein by reference. A specific algorithm for performing the basic water- 
pouring process for a MIMO-OFDM system is described in U.S. Patent Application 
Serial No. 09/978,337, entitled "Method and Apparatus for Determining Power 
Allocation in a MIMO Communication System," filed October 15, 2001. Water- 
pouring is also described in further detail in U.S. Patent Application Serial No. 
10/056,275, entitled "Reallocation of Excess Power for Full Channel-State Information 
(CSI) Multiple-Input, Multiple-Output (MIMO) Systems," filed January 23, 2002. 
These applications are assigned to the assignee of the present application and 
incorporated herein by reference. 

[1085] If water-pouring is performed to allocate the total available transmit power 
to the eigenmodes, then water-pouring analysis unit 234 provides a set of N s scaling 
values, b = {b 0 ... b, ... b Ns ], for the N s eigenmodes. Each scaling value is for a 
respective eigenmode and is used to scale the set of weights determined for that 
eigenmode. 

[1086] For eigenmode i, a set of weights, w l7 = [w..(0) ... w„(k) ... w ti (N F -l)] r , 
used to invert the channel and scale the transmit power of the eigenmode may be 
derived as: 
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= -f&- , for Jk = (0, 1, ... , N F -1) , Eq (10) 

A u (k) 

where the normalization factor, a n and the scaling value, b i , are derived as described 
above. 

[1087] Weight computation unit 230 provides the set of weighting matrices, W, 
which may be obtained using the weights w. (k) or w u (&) . Each weighting matrix, 
W(ifc), is a diagonal matrix whose diagonal elements are composed of the weights 
derived above. In particular, if only channel inversion is performed, then each 
weighting matrix, W(&) , for k = (0, 1, ... , N F -1) , is defined as: 

W(*) = diag (w n (k), w 22 (k), ... ,w NsNs (k)) , Eq (11a) 

where w u (k) is derived as shown in equation (3). And if both channel inversion and 
water-pouring are performed, then each weighting matrix, W(fc), for 
k = (0, 1, ... , N F — 1) , is defined as: 

W (*) = diag (*„(*), *„(*), ... ,*„,„,(*)) , Eq(llb) 

where w„(k) is derived as shown in equation (10). 

[1088] Referring back to FIG. 2, a scaler/IFFT 236 receives (1) the set of unitary 
matrices, V , which are the matrices of right eigen-vectors of H , and (2) the set of 
weighting matrices, W, for all N F frequency bins. Scaler/IFFT 236 then derives a 
spatio-temporal pulse-shaping matrix, P w («) , for the transmitter based on the received 
matrices. Initially, the square root of each weighting matrix, W(&), is computed to 
obtain a corresponding matrix, -Jwjk) , whose elements are the square roots of the 
elements of W(&). The elements of the weighting matrices, W(fc) for 
k = (0, 1, ... , N F -1) , are related to the power of the eigenmodes. The square root then 
transforms the power to the equivalent signal scaling. For each frequency bin k, the 
product of the square-root weighting matrix, <sJW(k) , and the corresponding unitary 
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matrix, V(fc) , is then computed to provide a product matrix, V(k)^JW(kj . The set of 
product matrices, \(k)JW(Jc) for k = (0, 1, ... ,N F -1), which is also denoted as 
V^W , defines the optimal or near-optimal spatio-spectral shaping to be applied to the 
modulation symbol vectors, s(n) . 

[1089] An inverse FFT of V^/W is then computed to derive the spatio-temporal 
pulse-shaping matrix, P lx (n) , for the transmitter, which may be expressed as: 

P lx (n) = WFT[VjW] . Eq(12) 

The pulse-shaping matrix, P lx (n) , is an N T xN T matrix. Each element of J* lx (n) is a 
set of Np time-domain values, which is obtained by the inverse FFT of a set of values 
for the corresponding element of the product matrices, V^W . Each column of J* lx (n) 
is a steering vector for a corresponding element of s(n) . 

[1090] A convolver 240 receives and preconditions the modulation symbol vectors, 
s(n) , with the pulse-shaping matrix, P„(n) , to provide the transmitted symbol vectors, 
x(n) . In the time domain, the preconditioning is a convolution operation, and the 
convolution of s(n) with P te (n) may be expressed as: 

x(") = 2XWs(n-^) • Eq(13) 

The matrix convolution shown in equation (13) may be performed as follows. To 
derive the i-th element of the vector x(n) for time n, x, (n) , the inner product of the i-th 
row of the matrix P„(^) with the vector s(n-£) is formed for a number of delay 
indices (e.g., 0<£< (N F -1)), and the results are accumulated to derive the element 
jc, («). The preconditioned symbol streams transmitted on each transmit antenna (i.e., 
each element of x(n) or *,.(«)) is thus formed as a weighted combination of the Nr 
modulation symbol streams, with the weighting determined by the appropriate column 
of the matrix P tx (n) . The process is repeated such that each element of x(n) is 
obtained from a respective column of the matrix P,,(") and the vector s(n) . 



WO 2004/002047 



PCT7US2003/019464 



24 

[1091] Each element of x(n) corresponds to a sequence of preconditioned symbols 
to be transmitted over a respective transmit antenna. The N T preconditioned symbol 
sequences collectively form a sequence of vectors, which are also referred to as the 
transmitted symbol vectors, x(n), with each such vector including up to N T 
preconditioned symbols to be transmitted from up to N T transmit antennas for the n-th 
symbol period. The N T preconditioned symbol sequences are provided to transmitters 
122a through 122t and processed to derive N T modulated signals, which are then 
transmitted from antennas 124a through 124t, respectively. 

[1092] The embodiment shown in FIG. 2 performs time-domain beam-steering of 
the modulation symbol vectors, s(n) . The beam-steering may also be performed in the 
frequency domain. This can be done using techniques, such as the overlap-add method, 
which are well-known in the digital signal processing field, for implementing finite- 
duration impulse response (FIR) filters in the frequency domain. In this case, the 
sequences that make up the elements of the matrix P tt (n) for n = (0, 1, ... ,N F -l) 
are each padded with N 0 - N F zeros, resulting in a matrix of zero-padded sequences, 
q^O*), for n = (0, 1, ... ,N 0 -l). An N Q -point fast Fourier transform (FFT) is then 
computed for each zero-padded sequence in the matrix q^C"), resulting in a matrix 
Q w (*) for * = ((), 1, ... ,N 0 -l). 

[1093] Also, the sequences of modulation symbols that make up the elements of 
s(n) are each broken up into subsequences of length N ss = N Q - N F + 1 . Each 
subsequence is then padded with N F -l zeros to provide a corresponding vector of 
length N 0 . The sequences of s(n) are thus processed to provide sequences of vectors 
of length N 0 , s e (n) , where the subscript £ is the index for the vectors that correspond 
to the zero-padded subsequences. An N 0 -point fast Fourier transform is then computed 
for each of the zero-padded subsequences, resulting in a sequence of frequency-domain 
vectors, S e (k), for different values of £. Each vector S e (k), for a given £, includes a 
set of frequency-domain vectors of length N 0 (i.e., for k = (Q, 1, ... ,N 0 -l)). The 
matrix Q tx (k) is then multiplied with the vector S e (k) , for each value of £ , where the 
pre-multiplication is performed for each value of k, i.e., for k = (0, 1, ... ,N 0 -l). 
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The inverse FFTs are then computed for the matrix-vector product Q tx (k)S e (k) to 
provide a set of time-domain subsequences of length N 0 . The resulting subsequences 
are then reassembled, according to the overlap-add method, or similar means, as is well- 
known in the art, to form the desired output sequences. 

[1094] FIG. 6 is a flow diagram of an embodiment of a process 600 that may be 
performed at the transmitter unit to implement the various transmit processing 
techniques described herein. Initially, data to be transmitted (i.e., the information bits) 
is processed in accordance with a particular processing scheme to provide a number of 
streams of modulation symbols (step 612). As noted above, the processing scheme may 
include one or more coding schemes and one or more modulation schemes (e.g., a 
separate coding and modulation scheme for each modulation symbol stream). 
[1095] An estimated channel response matrix for the MTMO channel is then 
obtained (step 614). This matrix may be the estimated channel impulse response 
matrix, £[, or the estimated channel frequency response matrix, H, which may be 
provided to the transmitter from the receiver. The estimated channel response matrix is 
then decomposed (e.g., using channel eigen-decomposition) to obtain a set of matrices 
of right eigen- vectors, V , and a set of matrices of singular values, A (step 616). 
[1096] A number of sets of weights, w l( , are then derived based on the matrices of 
singular values (step 618). One set of weight may be derived for each eigenmode used 
for data transmission. These weights are used to reduce or minimize intersymbol 
interference at the receiver by inverting the frequency response of each eigenmode 
selected for use. 

[1097] A set of scaling values, b, may also be derived based on the matrices of 
singular values (step 620). Step 620 is optional, as indicated by the dashed box for step 
620 in FIG. 6. The scaling values may be derived using water-pouring analysis and are 
used to adjust the transmit powers of the selected eigenmodes. 

[1098] A pulse-shaping matrix, P„(«), is then derived based on the matrices of 
right eigen-vectors, V, the sets of weights, w fI ., and (if available) the set of scaling 
values, b (step 622). The streams of modulation symbols are then preconditioned (in 
either the time domain or frequency domain) based on the pulse-shaping matrix to 
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provide a number of streams of preconditioned symbols, x(n) , to be transmitted over 
the MEVIO channel (step 624). 

[1099] Time-domain transmit processing with channel eigenmode decomposition 
and water-pouring is described in further detail in U.S. Patent Application Serial No. 
10/017,038, entitled 'Time-Domain Transmit and Receive Processing with Channel 
Eigen-mode Decomposition for MIMO Systems," filed December 7, 2001, which is 
assigned to the assignee of the present application and incorporated herein by reference. 
[1100] FIG. 7 is a block diagram of an embodiment of a receiver unit 700 capable 
of implementing various processing techniques described herein. Receiver unit 700 is 
an embodiment of the receiver portion of receiver system 150 in FIG. 1. Receiver unit 
700 includes (1) a RX MEVIO processor 160a that processes N R received symbol 
streams to provide N T recovered symbol streams and (2) a RX data processor 162a that 
demodulates, deinterleaves, and decodes the recovered symbols to provide decoded bits. 
RX MEVIO processor 160a and RX data processor 162a are one embodiment of RX 
MEMO processor 160 and RX data processor 162, respectively, in FIG. 1. 
[1101] Referring back to FIG. 1, the transmitted signals from Nt transmit antennas 
are received by each of Nr antennas 152a through 152r. The received signal from each 
antenna is routed to a respective receiver 154, which is also referred to as a front-end 
processor. Each receiver 154 conditions (e.g., filters, amplifies, and frequency 
downconverts) a respective received signal, and further digitizes the conditioned signal 
to provide ADC samples. Each receiver 154 may further data demodulate the ADC 
samples with a recovered pilot to provide a respective stream of received symbols. 
Receivers 154a through 154r thus provide N R received symbol streams. These streams 
collectively form a sequence of vectors, which are also referred to as the received 
symbol vectors, r(n) , with each such vector including Nr received symbols from the Nr 
receivers 154 for the n-th symbol period. The received symbol vectors, r(n) , are then 
provided to RX MEMO processor 160a. 

[1102] Within RX MEMO processor 160a, a channel estimator 712 receives the 
vectors r(n) and derives an estimated channel impulse response matrix, A , which may 
be sent back to the transmitter system and used in the transmit processing. An FFT 714 
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then performs an FFT on the estimated channel impulse response matrix, A , to obtain 
the estimated channel frequency response matrix, H (i.e., H = FFT [ti] ). 
[1103] A unit 716 then performs the channel eigen-decomposition of H(fc), for 
each frequency bin k, to obtain the corresponding matrix of left eigen-vectors, U(&) . 
Each column of U , where U = [U(0) ... U(fc) ... V(N F -1)] , is a steering vector for a 
corresponding element of r(n), and is used to orthogonalize the received symbol 
streams. An IFFT 718 then performs the inverse FFT of U to obtain a spatio-temporal 
pulse-shaping matrix, Uin) , for the receiver system. 

[1104] A convolver 720 then conditions the received symbol vectors, r(n) , with the 
conjugate transpose of the spatio-temporal pulse-shaping matrix, u" (n), to obtain the 
recovered symbol vectors, s(n) , which are estimates of the modulation symbol vectors, 
s(n) . In the time domain, the conditioning is a convolution operation, which may be 
expressed as: 

s(») = X^"^ (n ~ £) • E( 5( 14 > 

[1105] The pulse-shaping at the receiver may also be performed in the frequency 
domain, similar to that described above for the transmitter. In this case, the N R 
sequences of received symbols for the N R receive antennas, which make up the sequence 
of received symbol vectors, r(n) , are each broken up into subsequences of Nss received 
symbols, and each subsequence is zero-padded to provide a corresponding vector of 
length N 0 . The N R sequences of r(n) are thus processed to provide N R sequences of 
vectors of length N 0 , r,(n), where the subscript I is the index for the vectors that 
correspond to the zero-padded subsequences. Each zero-padded subsequence is then 
transformed via an FFT, resulting in a sequence of frequency-domain vectors, R e (k) , 
for different values of I. Each vector R e (k), for a given i, includes a set of 
frequency-domain vectors of length N Q (i.e., for k = (0, 1, ... ,N 0 -\) ). 
[1106] The conjugate transpose of the spatio-temporal pulse-shaping matrix, 
u" (n) , is also zero-padded and transformed via an FFT to obtain a frequency-domain 
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matrix, V" (k) for k = (0, 1, ... , N D -1) . The vector R e (k) , for each value of I , is 
then pre-multiplied with the conjugate transpose matrix ff (k) (where the pre- 
multiplication is performed for each value of k, i.e., for k = (0, 1, ... ,N 0 -l)) to 
obtain a corresponding frequency-domain vector S e (k). Each vector S e (k), which 
includes a set of frequency-domain vectors of length N 0 , can then be transformed via 
an inverse FFT to provide a corresponding set of time-domain subsequences of length 
N 0 . The resulting subsequences are then reassembled according to the overlap-add 
method, or similar means, as is well-known in the art, to obtain sequences of recovered 
symbols, which corresponds to the set of recovered symbol vectors, s(n) . 
[1107] Thus recovered symbol vectors, s(n) , may be characterized as a convolution 
in the time domain, as follows: 

s(n) = £r(£)s(n-^) + z(n) , Eq (15) 

where T(£) is the inverse FFT of A(k) = A(ifc)VW(fc) ; and 

z(«) is the received noise as transformed by the receiver spatio-temporal pulse- 
shaping matrix, u" (£) . 
The matrix Y(n) is a diagonal matrix of eigen-pulses derived from the set of matrices 
A, where A = [A(0) ... A(k) ... A(iV F -1)]. In particular, each diagonal element of 
T(n) corresponds to an eigen-pulse that is obtained as the IFFT of a set of singular 
values, [i,,(0) ... X u {k) ... X U {N F - 1)] T , for a corresponding element of A, 
[1108] The two forms for ordering the singular values, sorted and random-ordered, 
result in two different types of eigen-pulses. For the sorted form, the resulting eigen- 
pulse matrix, T s (n) , is a diagonal matrix of pulses that are sorted in descending order of 
energy content. The pulse corresponding to the first diagonal element of the eigen-pulse 
matrix, {r j (n)} n , has the most energy, and the pulses corresponding to elements further 
down the diagonal have successively less energy. Furthermore, when the SNR is low 
enough that water-pouring results in some of the frequency bins having little or no 
energy, the energy is taken out of the smallest eigen-pulses first. Thus, at low SNRs, 
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one or more of the eigen-pulses may have little or no energy. This has the advantage 
that at low SNRs, the coding and modulation are simplified through the reduction in the 
number of orthogonal subchannels. However, in order to approach channel capacity, 
the coding and modulation are performed separately for each eigen-pulse. 
[1109] The random-ordered form of the singular values in the frequency domain 
may be used to further simplify coding and modulation (i.e., to avoid the complexity of 
separate coding and modulation for each element of the eigen-pulse matrix). In the 
random-ordered form, for each frequency bin, the ordering of the singular values is 
random instead of being based on their magnitude or size. This random ordering can 
result in approximately equal energy in all of the eigen-pulses. When the SNR is low 
enough to result in frequency bins with little or no energy, these bins are spread 
approximately evenly among the eigenmodes so that the number of eigen-pulses with 
non-zero energy is the same independent of SNR. At high SNRs, the random-order 
form has the advantage that all of the eigen-pulses have approximately equal energy, in 
which case separate coding and modulation for different eigenmodes are not required. 
[1110] If the response of the MIMO channel is frequency selective, then the 
elements in the diagonal matrices, A(k) , are time-dispersive. However, because of the 
pre-processing at the transmitter to invert the channel, the resulting recovered symbol 
sequences, |(n), have little intersymbol interference, if the channel inversion is 
effectively performed. In that case, additional equalization would not be required at the 
receiver to achieve high performance. 

[1111] If the channel inversion is not effective (e.g., due to an inaccurate estimated 
channel frequency response matrix, H) then an equalizer may be used to equalize the 
recovered symbols, s(n) , prior to the demodulation and decoding. Various types of 
equalizer may be used to equalize the recovered symbol streams, including a minimum 
mean square error linear equalizer (MMSE-LE), a decision feedback equalizer (DFE), a 
maximum likelihood sequence estimator (MLSE), and so on. 

[1112] Since the orthogonalization process at the transmitter and receiver results in 
decoupled (i.e., orthogonal) recovered symbol streams at the receiver, the complexity of 
the equalization required for the decoupled symbol streams is greatly reduced. In 
particular, the equalization may be achieved by parallel time-domain equalization of the 
independent symbol streams. The equalization may be performed as described in the 
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aforementioned U.S. Patent Application Serial No. 10/017,038, and in U.S. Patent 
Application Serial No. 09/993,087, entitled "Multiple-Access Multiple-Input Multiple- 
Output (MIMO) Communication System," filed November 6, 2001, which is assigned 
to the assignee of the present application and incorporated herein by reference. 
[1113] For the embodiment in FIG. 7, the recovered symbol vectors, s(n) , are 
provided to RX data processor 162a. Within processor 162a, a symbol demapping 
element 732 demodulates each recovered symbol in in accordance with a 

demodulation scheme that is complementary to the modulation scheme used for that 
symbol at the transmitter system. The demodulated data from symbol demapping 
element 732 is then de-interleaved by a deinterleaver 734. The deinterleaved data is 
further decoded by a decoder 736 to obtain the decoded bits, d t , which are estimates of 
the transmitted information bits, d i . The deinterleaving and decoding are performed in 
a manner complementary to the interleaving and encoding, respectively, performed at 
the transmitter system. For example, a Turbo decoder or a Viterbi decoder may be used 
for decoder 736 if Turbo or convolutional coding, respectively, is performed at the 
transmitter system. 

[1114] FIG. 8 is a flow diagram of a process 800 that may be performed at the 
receiver unit to implement the various receive processing techniques described herein. 
Initially, an estimated channel response matrix for the MIMO channel is obtained (step 
812). This matrix may be the estimated channel impulse response matrix, fi, or the 
estimated channel frequency response matrix, H. The matrix A or H may be 
obtained, for example, based on pilot symbols transmitted over the MIMO channel. The 
estimated channel response matrix is then decomposed (e.g., using channel eigen- 
decomposition) to obtain a set of matrices of left eigen-vectors, U (step 814). 
[1115] A pulse-shaping matrix U(n) is then derived based on the matrices of left 
eigen-vectors, U (step 816). The streams of received symbols are then conditioned (in 
either the time domain or frequency domain) based on the pulse-shaping matrix U(n) 
to provide the streams of recovered symbols (step 818). The recovered symbols are 
further processed in accordance with a particular receive processing scheme, which is 
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complementary to the transmit processing scheme used at the transmitter, to provide the 
decoded data (step 820). 

[1116] Time-domain receive processing with channel eigenmode decomposition is 
described in further detail in the aforementioned U.S. Patent Application Serial No. 
10/017,038. 

[1117] The techniques for processing a data transmission at a transmitter and a 
receiver described herein may be implemented in various wireless communication 
systems, including but not limited to MIMO and CDMA systems. These techniques 
may also be used for the forward link and/or the reverse link. 

[1118] The techniques described herein to process a data transmission at the 
transmitter and receiver may be implemented by various means. For example, these 
techniques may be implemented in hardware, software, or a combination thereof. For a 
hardware implementation, the elements used to perform various signal processing steps 
at the transmitter (e.g., to code and modulate the data, decompose the channel response 
matrix, derive the weights to invert the channel, derive the scaling values for power 
allocation, derive the transmitter pulse-shaping matrix, precondition the modulation 
symbols, and so on) or at the receiver (e.g., to decompose the channel response matrix, 
derive the receiver pulse-shaping matrix, condition the received symbols, demodulate 
and decode the recovered symbols, and so on) may be implemented within one or more 
application specific integrated circuits (ASICs), digital signal processors (DSPs), digital 
signal processing devices (DSPDs), programmable logic devices (PLDs), field 
programmable gate arrays (FPGAs), processors, controllers, micro-controllers, 
microprocessors, other electronic units designed to perform the functions described 
herein, or a combination thereof. 

[1119] For a software implementation, some or all of the signal processing steps at 
each of the transmitter and receiver may be implemented with modules (e.g., 
procedures, functions, and so on) that perform the functions described herein. The 
software codes may be stored in a memory unit (e.g., memories 132 and 172 in FIG. 1) 
and executed by a processor (e.g., controllers 130 and 170). The memory unit may be 
implemented within the processor or external to the processor, in which case it can be 
communicatively coupled to the processor via various means as is known in the art. 
[1120] The previous description of the disclosed embodiments is provided to enable 
any person skilled in the art to make or use the present invention. Various 
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modifications to these embodiments will be readily apparent to those skilled in the art, 
and the generic principles defined herein may be applied to other embodiments without 
departing from the spirit or scope of the invention. Thus, the present invention is not 
intended to be limited to the embodiments shown herein but is to be accorded the widest 
scope consistent with the principles and novel features disclosed herein. 



[1121] WHAT IS CLAIMED IS: 
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CLAIMS 

1. In a multiple-input multiple-output (MIMO) communication system, a 
method for processing data for transmission over a MIMO channel, comprising: 

processing data in accordance with a particular processing scheme to provide a 
plurality of streams of modulation symbols; 

deriving a pulse-shaping matrix based on an estimated response of the MIMO 
channel and in a manner to reduce intersymbol interference at a receiver; and 

preconditioning the plurality of modulation symbol streams based on the pulse- 
shaping matrix to provide a plurality of streams of preconditioned symbols for 
transmission over the MIMO channel. 

2. The method of claim 1, further comprising: 

deriving a plurality of weights based on an estimated channel response matrix 
for the MIMO channel, wherein the weights are used to invert a frequency response of 
the MIMO channel, and wherein the pulse-shaping matrix is further derived based on 
the weights. 

3. The method of claim 2, further comprising: 

decomposing the estimated channel response matrix to obtain a plurality of 
matrices of ei gen- vectors and a plurality of matrices of singular values, and 

wherein the weights are derived based on the matrices of singular values and the 
pulse-shaping matrix is further derived based on the matrices of eigen-vectors. 

4. The method of claim 2, wherein the estimated channel response matrix is 
descriptive of a plurality of eigenmodes of the MIMO channel. 

5. The method of claim 4, wherein one set of weights is derived for each 
eigenmode used for data transmission and wherein the weights in each set are derived to 
invert the frequency response of the corresponding eigenmode. 

6. The method of claim 4, further comprising: 
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deriving a plurality of scaling values based on the matrices of singular values, 
wherein the scaling values are used to adjust transmit powers for the eigenmodes of the 
MEVIO channel, and wherein the pulse-shaping matrix is further derived based on the 
scaling values. 

7. The method of claim 6, wherein the scaling values are derived based on 
water-pouring analysis. 

8. The method of claim 3, wherein the estimated channel response matrix is 
provided in the frequency domain and is decomposed in the frequency domain. 

9. The method of claim 3, wherein the estimated channel response matrix is 
decomposed using channel eigen-decomposition. 

10. The method of claim 4, wherein eigenmodes associated with 
transmission capabilities below a particular threshold are not used for data transmission. 

1 1 . The method of claim 3, wherein the singular values in each matrix are 
sorted based on their magnitude. 

12. The method of claim 4, wherein the singular values in each matrix are 
randomly ordered such that the eigenmodes of the MLMO channel are associated with 
approximately equal transmission capabilities. 

13. The method of claim 1, wherein the pulse-shaping matrix comprises a 
plurality of sequences of time-domain values, and wherein the preconditioning is 
performed in the time domain by convolving the streams of modulation symbols with 
the pulse-shaping matrix. 

14. The method of claim 1, wherein the pulse-shaping matrix comprises a 
plurality of sequences of frequency-domain values, and wherein the preconditioning is 
performed in the frequency domain by multiplying a plurality of streams of transformed 
modulation symbols with the pulse-shaping matrix. 
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15. The method of claim 1, wherein the pulse-shaping matrix is derived to 
maximize capacity by allocating more transmit power to eigenmodes of the MIMO 
channel having greater transmission capabilities. 

16. The method of claim 1, wherein the pulse-shaping matrix is derived to 
provide approximately equal received signal-to-noise-and-interference ratios (SNRs) for 
the plurality of modulation symbol streams at the receiver. 

17. The method of claim 1, wherein the particular processing scheme defines 
a separate coding and modulation scheme for each modulation symbol stream. 

18. The method of claim 1, wherein the particular processing scheme defines 
a common coding and modulation scheme for all modulation symbol streams. 

19. In a multiple-input multiple-output (MIMO) communication system, a 
method for processing data for transmission over a MMO channel, comprising: 

processing data in accordance with a particular processing scheme to provide a 
plurality of streams of modulation symbols; 

obtaining an estimated channel response matrix for the MIMO channel; 

decomposing the estimated channel response matrix to obtain a plurality of 
matrices of eigen-vectors and a plurality of matrices of singular values; 

deriving a plurality of weights based on the matrices of singular values, wherein 
the weights are used to invert the frequency response of the MIMO channel; 

deriving a pulse-shaping matrix based on the matrices of eigen-vectors and the 
weights; and 

preconditioning the plurality of streams of modulation symbols based on the 
pulse-shaping matrix to provide a plurality of streams of preconditioned symbols for 
transmission over the MIMO channel. 

20. The method of claim 19, further comprising: 

deriving a plurality of scaling values based on the matrices of singular values, 
wherein the scaling values are used to adjust transmit powers for eigenmodes of the 
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MIMO channel, and wherein the pulse-shaping matrix is further derived based on the 
scaling values. 

21. A memory communicatively coupled to a digital signal processing 
device (DSPD) capable of interpreting digital information to: 

process data in accordance with a particular processing scheme to provide a 
plurality of streams of modulation symbols; 

derive a pulse-shaping matrix based on an estimated response of the MIMO 
channel and in a manner to reduce intersymbol interference at a receiver; and 

precondition the plurality of streams of modulation symbols based on the pulse- 
shaping matrix to provide a plurality of streams of preconditioned symbols for 
transmission over the MIMO channel. 

22. In a multiple-input multiple-output (MIMO) communication system, a 
method for processing a data transmission received via a MIMO channel, comprising: 

obtaining an estimated channel response matrix for the MIMO channel; 

decomposing the estimated channel response matrix to obtain a plurality of 
matrices of eigen-vectors; 

deriving a pulse-shaping matrix based on the matrices of eigen-vectors; and 

conditioning a plurality of streams of received symbols based on the pulse- 
shaping matrix to provide a plurality of streams of recovered symbols which are 
estimates of modulation symbols transmitted for the data transmission, wherein the 
modulation symbols are preconditioned at a transmitter, prior to transmission over the 
MIMO channel, in a manner to reduce intersymbol interference at a receiver. 

23. The method of claim 22, wherein the conditioning is performed in the 
time domain based on a time-domain pulse-shaping matrix. 

24. The method of claim 22, wherein the conditioning is performed in the 
frequency domain and includes 

transforming the plurality of received symbol streams to the frequency domain; 
multiplying the transformed received symbol streams with a frequency-domain 
pulse-shaping matrix to provide a plurality of conditioned symbol streams; and 
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transforming the plurality of conditioned symbol streams to the time domain to 
provide the plurality of recovered symbol streams. 

25. The method of claim 22, wherein the conditioning orthogonalizes a 
plurality of streams of modulation symbols transmitted over the MTMO channel. 

26. The method of claim 22, further comprising: 

demodulating the plurality of recovered symbol streams in accordance with one 
or more demodulation schemes to provide a plurality of demodulated data streams; and 

decoding the plurality of demodulated data streams in accordance with one or 
more decoding schemes to provide decoded data. 

27. The method of claim 22, further comprising: 

deriving channel state information (CSI) comprised of the estimated channel 
response matrix for the MTMO channel; and 

sending the CSI back to the transmitter. 

28. In a multiple-input multiple-output (MEVIO) communication system, a 
method for processing a data transmission received via a MIMO channel, comprising: 

obtaining an estimated channel response matrix for the MTMO channel; 

decomposing the estimated channel response matrix to obtain a plurality of 
matrices of ei gen- vectors; 

deriving a pulse-shaping matrix based on the matrices of eigen-vectors; 

conditioning a plurality of streams of received symbols based on the pulse- 
shaping matrix to provide a plurality of streams of recovered symbols which are 
estimates of modulation symbols transmitted for the data transmission, wherein the 
modulation symbols are preconditioned at a transmitter, prior to transmission over the 
MTMO channel, in a manner to reduce intersymbol interference at a receiver; 

demodulating the plurality of recovered symbol streams in accordance with one 
or more demodulation schemes to provide a plurality of demodulated data streams; and 

decoding the plurality of demodulated data streams in accordance with one or 
more decoding schemes to provide decoded data. 
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29. A memory communicatively coupled to a digital signal processing 
device (DSPD) capable of interpreting digital information to: 

obtain an estimated channel response matrix for a MTMO channel used for a data 
transmission; 

decompose the estimated channel response matrix to obtain a plurality of 
matrices of ei gen- vectors; 

derive a pulse-shaping matrix based on the matrices of eigen-vectors; and 
condition a plurality of streams of received symbols based on the pulse-shaping 
matrix to provide a plurality of streams of recovered symbols which are estimates of 
modulation symbols transmitted for the data transmission, wherein the modulation 
symbols are preconditioned at a transmitter, prior to transmission over the MIMO 
channel, in a manner to reduce intersymbol interference at a receiver. 

30. A transmitter unit in a multiple-input multiple-output (MIMO) 
communication system, comprising: 

a TX data processor operative to process data in accordance with a particular 
processing scheme to provide a plurality of streams of modulation symbols; and 

a TX MIMO processor operative to derive a pulse-shaping matrix based on an 
estimated response of a MTMO channel and in a manner to reduce intersymbol 
interference at a receiver, and to precondition the plurality of modulation symbol 
streams based on the pulse-shaping matrix to provide a plurality of streams of 
preconditioned symbols for transmission over the MIMO channel. 

3 1 . The transmitter unit of claim 30, wherein the TX MIMO processor is 
further operative to derive a plurality of weights based on an estimated channel response 
matrix for the MIMO channel, wherein the weights are used to invert the frequency 
response of the MIMO channel, and wherein the pulse-shaping matrix is derived based 
in part on the weights. 

32. The transmitter unit of claim 3 1 , wherein the TX MIMO processor is 
further operative to decompose the estimated channel response matrix to obtain a 
plurality of matrices of eigen-vectors and a plurality of matrices of singular values, and 
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wherein the weights are derived based on the matrices of singular values and the pulse- 
shaping matrix is further derived based on the matrices of eigen-vectors. 

33. The transmitter unit of claim 3 1 , wherein the TX MTMO processor is 
further operative to derive a plurality of scaling values used to adjust transmit powers 
for the eigenmodes of the MTMO channel, and wherein the pulse-shaping matrix is 
further derived based on the scaling values. 

34. The transmitter unit of claim 33, wherein the scaling values are derived 
based on water-pouring analysis on a plurality of matrices of singular values obtained 
from the estimated channel response matrix. 

35. An apparatus in a multiple-input multiple-output (MEMO) 
communication system, comprising: 

means for processing data in accordance with a particular processing scheme to 
provide a plurality of streams of modulation symbols; 

means for deriving a pulse-shaping matrix based on an estimated response of a 
MIMO channel and in a manner to reduce intersymbol interference at a receiver; and 

means for preconditioning the plurality of modulation symbol streams based on 
the pulse-shaping matrix to provide a plurality of streams of preconditioned symbols for 
transmission over the MIMO channel. 

36. A digital signal processor comprising: 

means for processing data in accordance with a particular processing scheme to 
provide a plurality of streams of modulation symbols; 

means for deriving a pulse-shaping matrix based on an estimated response of a 
multiple-input multiple-output (MIMO) channel and in a manner to reduce intersymbol 
interference at a receiver; and 

means for preconditioning the plurality of modulation symbol streams based on 
the pulse-shaping matrix to provide a plurality of streams of preconditioned symbols for 
transmission over the MIMO channel. 
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37. A receiver unit in a multiple-input multiple-output (MEVIO) 
communication system, comprising: 

an RX MMO processor operative to obtain an estimated channel response 
matrix for a MEVIO channel used for a data transmission, decompose the estimated 
channel response matrix to obtain a plurality of matrices of eigen- vectors, derive a 
pulse-shaping matrix based on the matrices of eigen-vectors, and condition a plurality of 
streams of received symbols based on the pulse-shaping matrix to obtain a plurality of 
streams of recovered symbols which are estimates of modulation symbols transmitted 
over the MIMO channel, wherein the modulation symbols were preconditioned at a 
transmitter, prior to transmission over the MEVIO channel, in a manner to reduce 
intersymbol interference at the receiver unit; and 

an RX data processor operative to process the plurality of recovered symbol 
streams in accordance with a particular processing scheme to provide decoded data. 

38. The receiver unit of claim 37, wherein the RX MIMO processor is 
operative to condition the plurality of streams of received symbols in the time domain 
based on a time-domain pulse-shaping matrix. 

39. An apparatus in a multiple-input multiple-output (MIMO) 
communication system, comprising: 

means for obtaining an estimated channel response matrix for a MEMO channel 
used for a data transmission; 

means for decomposing the estimated channel response matrix to obtain a 
plurality of matrices of eigen-vectors; 

means for deriving a pulse-shaping matrix based on the matrices of eigen- 
vectors; and 

means for conditioning a plurality of streams of received symbols based on the 
pulse-shaping matrix to provide a plurality of streams of recovered symbols which are 
estimates of modulation symbols transmitted for the data transmission, wherein the 
modulation symbols are preconditioned at a transmitter, prior to transmission over the 
MEVIO channel, in a manner to reduce intersymbol interference at a receiver. 



40. 



A digital signal processor comprising: 
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means for obtaining an estimated channel response matrix for a multiple-input 
multiple-output (MIMO) channel used for a data transmission; 

means for decomposing the estimated channel response matrix to obtain a 
plurality of matrices of eigen-vectors; 

means for deriving a pulse-shaping matrix based on the matrices of eigen- 
vectors; and 

means for conditioning a plurality of streams of received symbols based on the 
pulse-shaping matrix to provide a plurality of streams of recovered symbols which are 
estimates of modulation symbols transmitted for the data transmission, wherein the 
modulation symbols are preconditioned at a transmitter, prior to transmission over the 
MIMO channel, in a manner to reduce intersymbol interference at a receiver. 
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